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Abstract 

Possible field-temperature superconducting (SC) phase diagrams in two-band quasi-two- 
dimensional materials with a strong paramagnetic pair-breaking (PPB) are considered theoreti¬ 
cally. Attention is paid to the case under a magnetic field perpendicular to the SC layers and 
to essential differences from the counterpart in the ordinary single-band material. It is found by 
examining the Hc 2 {T) curve and the vortex lattices to be realized close to Hc 2 that a PPB-induced 
SC phase with a spatial modulation parallel to the field tends to occur more easily than in the 
single-band case, and that a crisscrossing vortex lattice proposed previously can occur in place of 
the conventional Fulde-Ferrell-Larkin-Ovchinnikov state within a parameter range. The relevance 
of the obtained results to FeSe and other iron-based superconductors is discussed. 
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1. Introduction 


Over the last decade, superconducting (SC) materials showing remarkably strong para¬ 
magnetic pair-breaking (PPB) effects have been discovered, and studying the PPB effect 
in superconductors with unconventional electronic properties is now one of the important 
theoretical subjects in condensed matter physics. In the quasi-two-dimensional (2D) heavy- 


fermion superconductors CeCoInj 


[l| and NpPd 2 Al 5 j^, the SC transition at high fields and 
hence at low temperatures has been shown to be nearly discontinuous, and, on the basis 
of this and the temperature dependence of the resulting Hc^iT) curve, an extremely strong 
’^PB effect in these materials has been clarified. Furthermore, the strange high-field phase 
3|, found in CeCoIns in a field parallel to the SC layers has been interpreted, by an 
approach based on a single-band electronic model, as a vortex lattice corresponding to the 
appearance of a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state ^ with a PPB-induced 


=pat.l .odulatioa paaalle. .o .he n.agae.ic field gQ. The appeaeaace of .fie PPB-a.daced 
spatial modulation parallel to the field rather than the corresponding perpendicular mod¬ 
ulation { 9 ] is a rare event in the weak-coupling BCS model and has been interpreted 
[n| as a reflection of the fact that this material at ambient pressure is close to an anti¬ 
ferromagnetic quantum critical point. Although the presence of a FFLO vortex state has 
also been suggested in CeCoIns in a field perpendicular to the a6-plane, i.e., the SC layers 
(H II c) [12|, one expects the FFLO state to usually appear in a field parallel to the SC 
layers in various quasi-2D materials including organics 13j. Several recent experiments on 
iron-based superconductors have also suggested SC properties induced by PPB in the fields 


parallel 


14 


Iq l and perpendicular [l6| to the SC layers such as a suppressed Hc 2 {T) curve 


at low temperatures. However, the iron-based superconductors are typically two-band ma¬ 
terials, and thus qualitatively different properties based on PPB may be expected to occur 
in these systems. In fact, clear evidence on the discontinuous ifc 2 transition has not been 
reported in iron-based materials, except for a few examples [l^. Besides this, the recent 
thermal conductivity measurement on the superconductor FeSe has led to the detection of 


a high-field SC phase in H 


Q. 


From the large Maki parameter estimated on one of two 


bands, it is reasonable to expect this high field SC phase to stem from a large PPB effect. 

In this work, we study possible field-temperature phase diagrams of two-band quasi-two- 
dimensional superconductors in a magnetic field perpendicular to the layers by constructing 
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the Ginzburg-Landau (GL) free-energy functional based on a microscopic model. In contrast 


to the single-band case 


10l | where the PPB strength is uniquely characterized by a single 


Maki parameter om, one cannot predict the resulting phase diagram only on the basis 
of the PPB strength in multi-band cases where a different ctm is dehned on a different 
band. By focusing on a two-band model in which the two oms remarkably differ from each 
other, we have a couple of categories of typical phase diagrams depending on the details 
or the breakdown of the parameters determining oms. When the interband coupling is 
substantial, we hnd two typical features that do not appear in the single-band case mentioned 
above. Firstly, we have a wide parameter range in which a PPB-induced spatial modulation 
characterizing a FFLO state is formed in the direction parallel to the held. In this state, the 
PPB-induced spatial modulation of a SG order parameter has the following form: A(r) oc 
exp(i/c^) for a FF state and A(r) oc cos{kz) for a LO state, where fc is a modulation 
wave number and the z-axis is parallel to the held. Secondly, over some range of material 
parameters, we hnd the situations in which the low-temperature Hc 2 transition is not a 
discontinuous but a continuous one to a vortex lattice with a modulation parallel to the held. 
This hnding seems to be consistent with the fact that, in the iron-based superconductors, 
the hrst-order-like ifc 2 transition has been rarely seen so far even in the cases indicative of a 
strong PPB ehect. On the basis of a previous argument given in relation to phenomena 
in GeGoIns, these two features satisfy the condition for the emergence of a crisscrossing 
vortex lattice [20[ as the PPB-induced modulated high-held SG state, where, in contrast to 
a FFLO state, the order parameter has only spatial line nodes. In addition, stable half- 
hux quanta can emerge in this state. Gonversely, if the interband coupling is sufficiently 
weak, the resulting phase diagram includes, as in the single-band case, a vortex state with a 
modulation in the plane perpendicular to the applied held at least within the weak-coupling 


BGS model 


io|. 


The outline of this manuscript is as follows. In Sect. 2, the model and methods of 
theoretical calculations are shown. In Sect. 3, the obtained phase diagrams are presented. In 
Sect. 4, our results are discussed in relation to the experimental observations. In Appendix, 
the details of derivations of theoretical expressions are explained. 
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2. Methods 


We start with the following quasi-2D two-band Hamiltonian "H in an external held H 
parallel to the c-axis: 


"H = 7^20 + 'Hhop + 'Hint 


Nu 


= Y. Y. E [<(r>^')]H^nHV-eA(r)]-aJn}C(r,j) 


i=i 
Ni. 


n=l,2a=±l 


J''layer f J ( 

+ Y Y E oi [C(r,J + l)]Vn(r,j) + [C(r,i)]Vn(r,J + 1) 

j=l n=l,2o-=±l ^ I 


-1 -^layer 

- - 

^ i=i 

The hrst line H 2 D and the second line Hhop in Eq- dl]) describe the in-plane and inter-plane 
motions of quasiparticles in a magnetic held, respectively, while the third term Hint denotes 
mutual interactions between quasiparticles. Furthermore, we use the following symbols: H 
the area of each layer; Wayer the number of layers; £n(p) (n = 1 and 2 ) the energy dispersion 
of a quasiparticle in the a 6 -plane on the n-th band; H = gn^JiBH/2 the Zeeman energy on the 
n-th band; a the spin projection index; gn the ehective ^f-factor on each band; /tb the Bohr 
magneton; J the interlayer hopping energy; and > 0 and (H)i 2 = (H) 2 i = M 3 the 

intraband and interband interaction constants. Note that the 2 ;-axis is parallel to the c-axis, 
and that r and j are the in-plane coordinate and layer number, respectively. Concerning 
each energy dispersion e:n(p), cylindrical symmetry is assumed for simplicity. The field 
operator of quasi-particles is dehned as V’n(r, j) = (HiViayer)”^'^^ Sp,fcexp[z(p-r + kj)]a'^ip,k) 
{—71 < k < 77 ) with the corresponding annihilation operator a‘^{p,k). Hereafter, the type 
H limit is assumed so that the magnetic held in the system is approximately equal to the 
applied held. At least near the resulting Hc 2 {T) line, which is the so-called upper critical 
held in the mean held approximation, this treatment is safely valid. In addition, except in 
the case with an extremely large Hc 2 { 0 ), we can use the quasi-classical approximation in 
which the ehect of the magnetic held on the one-particle Green’s function on the n-th band 
is represented in terms of a phase factor as 


dhJ2iU)nn' Y [C(r,j)]n^n"(r,j)]Vn"(r,i)C'(r,j)- 


( 1 ) 


^n(r,r', ~ exp 


ie J d\- A(l) 


1 

f^-ZVlayer 


Y(iW{i[p-{r-r') + k{j-j')])}g^{p,k,Ui), 

p,k 

( 2 ) 
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where ^^(p,/c,cjj) = [iui—^n{p)+o'In — Jcosk]~^,^n{p) = £n(p)—/W with chemical potential 
fi, and ui is the Fermion Matsubara frequency. Below, by the functional integral method 
and the Hubbard-Stratonovich (HS) transformation, we derive the GL functional and then 
examine the Hc 2 {T) line and the nature of the mean held transition in a two-band system. 
The ensuing treatment will be explained hereafter separately depending on the sign of |17| = 
U1U2 — 

(i) |G|>0 

When \U\ > 0, as shown in Appendix, the auxiliary held A„(r,j) corresponding to the 
SC order parameter held j)) can be straightforwardly introduced 

to obtain the GL functional 

Ggl [Ai, A 2 ] = G(o) (U-^) , A„, 

3 Ti^n' 



where r2(o) is the thermodynamic potential in the normal phase, 

S' = - drY, J d^Y An^I(r)^i(r) + A;^/>;^(r)^/>;^(r) , (4) 

and denotes the ensemble average of X with respect to the noninteracting part of 

the action with only the connected diagrams retained. We note that here the Fermion helds 
r) and r) are the Grassmann numbers. 

First, let us derive the quadratic ((T(|Ap)) terms. In this case, the quadratic term in 
Ggl takes the form 

^(2) = E E / d^T k) - 5nn'i?{2)n(n, k)] k), (5) 

n,n' k 

where 

= Y. Ys:{p,k',^,)5^’{n-p,k- (6) 

^“JViayer p,a;; o-=±l k' 

n = —iV —2eA(r) is the gauge-invariant gradient operator, and the Fourier transformation 
is dehned as A„(r,j) = (Wayer)”^^^ 6xp(*^j)An(i',/i^)- Following the route sketched in 
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Appendix, we introduce the parameter p and the cutoff pcut instead of the energy cutoff, 
and then we obtain 


n,n' A^=0 k 
oo „ 


-1 


^2) = EEE/^'r[AW(r,fc)]*|(t/ 

U2-D,f[^\k-T,H) 


5n,n'Dnfi^\k;T,H) 


k) 


= E E / 

N=0 k 


|AS^^(r, A ;)|2 + 


Ui-D,fi^\k-T,H) 


X |A^^V)^)I 


^ - Ms 


[AE(r, k) + (c. c.) 


W, 


where 


1*00 r 

fP\t,T,H) = 2nT dp- 

SI 


°° cos( 24 p) 


Pcut sinh(27rTp) 


k 

2Jpsin( - 


exp 


Vn^P^ 


4A 


H 


L 


N 


(vll_ 

\2Xh‘^ 


(7) 


, ( 8 ) 


A^^)(r, fc) represents the A-th Landau level component of A„(r, fc), Uj = Uj/\U\^ Jo(^) 
represents the zeroth-order Bessel function, and L 7 v(a^) represents the A-th-order Laguerre 
polynomial. Furthermore, Dn and are the density of states (DOS) at the Fermi surface 
and the Fermi velocity on the n-th band, respectively, and Xh = (2|e|A)“^/^ is the magnetic 
length. By diagonalizing the expression of Eq. ([7]) and using the resulting eigenvalues 
> cE) and the ‘eigenfunctions’ Ajf^^ = X]n ^( 2 ) may be expressed as 

00 . 

Ar=0 k 

We call A^^ merely as the ‘± component’ of Therefore, if the mean held transition is 

of the second order, the equation determining Hc 2 {T) is c^^\k-, T, H) = 0, or equivalently, 

1 - Ai/f)(fc;T, A) - A 2 /r(fc;T, A) + - As^)fi^\k;T, H)f^^\k;T, H) = 0, (10) 

where A and k should be chosen so that the magnitude of Hc 2 becomes highest for a given T, 
and Ai = DiUi, A 2 = D 2 U 2 , and A 3 = y/lA[D^u^ are the dimensionless coupling constants. 

Next, we turn to deriving the quartic ((!1(|A|^)) term in Dgl- Formally, this term is 
expressed as 

^(4) = fE E [ c?^r7F(4)„({ni,/ci})A;(ri,fci)A„(r2,fc2) 


A /Cl,/C 2 ,fc 3 

X A*(r3,/c3)A„(r4, fci + fcs - ^ 2 ) 


( 11 ) 


where 


■77(4)n ({nj, fcj}) 


T 


DA, 


layer p^uii 


E ((p> k, uji)g^^ (n* - p, fci - k, -ui) 


X ^„‘"(n2 - p,/C2 - -M;;)^)((n3 - 112 + P,/C3 - ^2 + fc,M;z)) ,(12) 


k^a 
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and {X)y denotes the average of X with respect to Y. Becanse we are interested in the 
region close to the resnlting Hc 2 {T), we only have to focns on the component’ dehned 
throngh the above-mentioned diagonalization, to examine the sign of the GL-qnartic term 
and thns the natnre of the mean held Hc 2 transition. Fnrthermore, for the moment, we will 
restrict ourselves to the lowest Landau level (LLL) modes of the order parameter, denoted 
as below, because, as will be explained in the next section, we often encounter the 
situations dominated by the LLL modes in the present two-band cases despite the presence 
of a large paramagnetic pair-breaking. We note that, in the familiar single-band case, a larger 
paramametic pair-breaking makes the vortex state a state dominated by higher Landau level 
modes [3, Q- By using the Landau gauge (A = Bxy), the order parameter will be written 
below in the form 


T, H) ^{o) ^ aL°^ (r, -k) 

\ / -A^lo-irpr 


A. 


layer 


W ^ n{k,T H) ^ ^ ^ 

layer S}J q 


(13) 


where Sh = y/^^uLy, Ly is the system 


size m 


the ^/-direction, and Uq{r) = 


exp 


— {x + qXn'^y /+ iqy ■ Using Eq. flT3l) and following the procedures to be ex¬ 


plained in Appendix, the quartic term is rewritten in the form 


Q 


(4) 


^ [w'^UM T, {<,.}, k; T, H) 4,-(qi,k) 


LLL, 


\/2Aiayer5'H nqgq 

X 0(^2, k)(j)*{q 3 , k)(j){q 4 , k) + K( {gj, fc; T, H) 0*(gi, k)(j){q 2 , k)(j)*{qs, -k)(j){q4, -k) 
+ K( k] T, H) 0*(gi, /c)0(g2, -k)(f)*{q^, -fc)0(g4, k) + {k ^ -k) 


(14) 


where = qi- qj, 


Vn = 27lTDr, 


COs2/„(pi +P2 +P3) . . , 3 . fk\ 

dpidp 2 dps . . . - ^^^Jo 2 J Pi + p 2 + P3 sm - 

smh 27rT(pi + P 2 + Pall V V2/ 


X (Reexp [5({gi,/3„,i})]|p4=o)^, 


14 = 2 'KTDn f dpidp 2 dp 3 Jo V + P 2 - Pa cos(fc))^ -F (pa sm{k)y 

Jo smh[27rT(pi-F p 2 psjj V 

X sinf^))(Reexp [R({gi,/3„,i})]|p4=o\ , (15) 


14 is Vn with Pi and po exchanged with each other, and B{{qi, /3n,i}) will be dehned in 
Appendix. If we assume that the k- and g*-dependences of both I 4 and W are weak and 
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negligible, the quartic term is expressed in the local form 


(~)local 

“(4) 


LLL,±fc 


= E|A<_“>(r,j)l*. ( 16 ) 


where 

= 27rTDr, 


' /Re 

^'sinh[27rT(pi+p2 + P3)] \ 




P4=0 


117) 


By using hl| 4 )'^^|LLL,±fc) we estimate the region in which a first-order transition occurs; if 
the quartic term is positive on the line where the quadratic term vanishes, a second-order 
transition occurs; otherwise, a hrst-order transition is expected to occur at a higher held. 


(ii) |C/|<0 


When \U\ <0, either of the two eigenvalues of U is negative, and then the interaction 
Hamiltonian takes the form 


n,n' 

where <h± are appropriate linear combinations of 'ipi'ipn- Because U- is negative, the second 
term of Eq. 0181) may be regarded as a contribution implying a repulsive interaction, and 
thus we neglect the term from now on. Then, as in the case of |f/| > 0, performing the 
HS transformation leads to a GL functional. Then, the corresponding quadratic term 17(2) 
[corresponding to Eq. ([9])] becomes 


n 


(2) - Z! Z! 

N=0 k 


U+ n 




(19) 


where w is the eigenvector that belongs to the eigenvalue and ■ We 

note that A^^ is different from A^^ dehned in the |f/| > 0 case. The quartic term 12(4) has 
the same form as in the |f/| > 0 case [see Eqs. flTTl) and flTTl) ]. if w^d]n{k]T,H) and A® in 
|t/| > 0 case are replaced by Wn and A+\ respectively. 


3. Results 

By numerically calculating the coefficients of the quadratic and quartic GL terms, we 
estimate the line in the case where the PPB effect is strong on one of the two bands. 
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Below, we will focus especially on this case by taking account of the SC properties seen in 
FeSe (see Sect. 1) 18j. We represent the relative strength of the paramagnetic effect to the 
orbital effect on each band via the parameter = CgnTco/ijriQVn^) oc Qn/vn^ proportional 
to the Maki parameter ctm, where C = 7 r/ 2 exp( 7 ) ~ 0.88, 7 is the Euler constant, ttiq is 
the electron mass, and T^o is the mean held transition temperature in the absence of the 
magnetic held. 

Let us start with reviewing the single-band limit in which A 2 = A 3 = 0, and the PPB 
on the only attractive band is given by ai = /(^tiTco)^ where is ^^£ 2 ( 0 ) in 

the single-band case with no PPB ehect included. In the weak PPB limit where ai <C 1, 
the orbital pair-breaking is dominant, and the second-order 11^2 transition to a conventional 
vortex lattice with no additional spatial modulation described in LLL occurs. Conversely, 
when the PPB is sufficiently strong with, say, ai > 1.4, the vortex lattice just below the 
Hc 2 {T) line is dominated by higher LLs at low but hnite temperatures at least within the 
weak-coupling BCS model although the Hc 2 transition is of the second order. This higher LL 
vortex state includes additional spatial modulations accompanied by antivortices in the plane 
perpendicular to the held. Depending on material parameters such as disorder strength, this 
higher LL transition is replaced by a discontinuous iLc 2 transition to the ordinary vortex 
lattice described in LLL, which may have a PPB-induced spatial modulation only in the 
direction parallel to the held. 

Hereafter, we will focus on the case in which PPB is strong on band 1 and weak on band 
2 by keeping the relation 02 = ai/16, and, regarding the intraband couplings, the values 
Ai = 0.2 and A 2 = 0.15 will be assumed. The interlayer hopping J will be hxed to lOvrTco 
in all of our computations. Furthermore, to reduce the number of material parameters 
appearing in our theoretical expressions, the in-plane Fermi momenta of the two bands are 
assumed to be the same so that D 2 /D 1 = where m* is the ehective mass on each 

band. Qualitatively, theoretical results are unahected by this simplihcation. Then, we have 
the relation a\loi 2 = giV 2 ^/{g 2 Vi^) = ginrdf'/{g 2 nfl^). Note that the ratio D 1 /D 2 can be 
changed while ai/a 2 is hxed. Hence, even if ai and 02 are hxed, the character of the Hc 2 
transition may be changed for diherent values of the DOSs Di and D 2 . 

Under such a parametrization, we have three control parameters that may ahect the 
resulting phase diagram, i.e., the interband coupling IA 3 I, the band asymmetry in the 

ehective mass (or the density of states), and another band asymmetry gi/g 2 in the Zeeman 
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Fig. 1. Hc 2 {T) (black solid) curve for m\ = gi = g2, and ai = 16a2 = 0.5. The coupling 
constants Ai = 0.2, A 2 = 0.15, and A 3 = 0.15 are used. The nature of the transition on Hc 2 {T) 
is of the second order at any field, and, because of the large asymmetry in the Fermi velocities on 
the two bands, the Hc 2 line has a positive curvature in intermediate fields. The A^ = 1 LL (black 
dotted) curve at which changes its sign is drawn for comparison and is not a transition line. 


term. Roughly speaking, (71/(72 determines the relative strength of only PPB between the 
two bands, while m^/m^ also affects the ratio of the coherence lengths, i.e., the orbital 
pair-breaking. 

First, we examine the case with a vanishingly weak interband coupling so that IA3I -C Ai 
and A2. Clearly, the SC transition in this limiting case is dominated by the states on one 
band. If the intraband couplings Ai and A2 are almost equal to each other, the band with 
a smaller PPB effect dominates in high fields. In the present case where the difference 
|Ai — A2I/A1 is of order unity, however, the situation is close to the above-mentioned single¬ 
band limit with A2 = A3 = 0, and the more attractive electronic states on band 1 largely 
determine the superconductivity despite the relation ai S> Q!2. In any case, the behaviors in 
the high-held and low-temperature region are similar to those in the single-band case 21 |. 
First, the transition character depends only on the value of ai and does not depend on its 
breakdown. Second, a hrst-order transition within the LLL or the transition to a vortex 
lattice dominated by a higher Landau level occurs at such low temperatures in the H-T 
phase diagram ff^ . 

In contrast, in the case with a moderately large interband coupling where IA3I ~ Ai or 
I A3 1 ~ A2, a two-band effect appears, and the character of the Hc2 transition depends not 
only on the value of an but also on their breakdown. First, we will explain the case in which 
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Fig. 2. (Color online) The Hc 2 {T) curve (black solid line with blue cross points) following from 
larger an values, ai = 16a2 = 3.0. Values of other parameters are the same as those in Fig. [U and 
thus the band asymmetry in arises from that of the effective mass. The Hc 2 {T) curve is defined 
in LLL. On its black line portion, c^^(O) changes its sign, and a second-order Hc 2 transition in the 
mean held approximation occurs between the conventional vortex lattice and the normal phase. 
On the other hand, the Hc 2 transition occurs as a hrst-order one in the region specihed by blue 
cross points. For reference, the N = 1 LL (black dotted) curve is also shown, as in Fig. [TJ 


the band asymmetry is due to that of the effective mass for the hxed interband coupling 
A 3 = 0.15. Note that this case belongs to the category with \U\ >0. In the limiting case with 
no PPB, the Hc 2 {T) curve in this case is well known and has an inflection in intermediate 


helds 


22l |. Such an inflection of Hc 2 {T) is also seen in Fig. [H where the weak PPB has been 


assumed by setting ai = 0.5. In this and the ensuing hgures, the temperature is measured 
in the unit of Tco, while the magnetic held is measured in the unit of IIc2(0). In each hgure, 
the black-solid line denotes the second-order transition line in the mean held approximation 
on which nucleation of the order parameter in LLL occurs, and the black-dotted curve is the 
corresponding one in the next lowest (N = 1) Landau level (LL) obtained by neglecting the 
presence of the SC ordering in LLL. In the present situation where the LLL ordering occurs 
at a higher held, no real transition occurs on this line associated with the N = 1 LL modes. 

Next, we show an example with larger PPB strengths in Fig. [2l As in the single-band 
case, the Hc 2 (T) curve takes a suppressed form. However, higher LL vortex states with a 
PPB-induced spatial modulation peroendicular to the held are not realized even in high helds 


in contrast to the single-band case 


aQ- 


Both the Hc 2 line and high-held vortex state are 


described in LLL. At the temperatures indicated by the blue cross points, the Hc 2 transition 
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is of the first order because the coefficient of the |A|^ term is negative there. Furthermore, 
the high-held vortex state might be accompanied by a spatial modulation parallel to the 
held. However, in Fig. m the region of a modulated phase has not been identihed. Further, 
we have verihed that such a phase diagram with a hrst-order Hc2 transition dehned within 
LLL is also obtained when IA3I is larger so that \U\ < 0 . It does not appear that the sign 
change of \U\ leads to a qualitative change in the resulting phase diagram. Although the Hc2 
transition under a similar situation where the interband coupling is dominant (IA3I 3 > Ai, A2) 
and the band asymmetry V2/V1 is large has been discussed elsewhere |22|, the possibility of 
the hrst-order Hc2 transition has not been examined there. 

Now, we will explain the case in which the band asymmetry is due primarily to that of 
the Zeeman energy. The dependences of the phase diagrams belonging to this case on the 
interband coupling IA3I are shown in Fig. [ 3 l The hgure (a) is a result for a relatively small 
A3 and, rehecting the result in the single-band case, the instability line to the N = 1 LL 
vortex state is close to the Hc2{T) curve dehned in LLL. As in Fig. | 2 l the Hc2 transition is 
discontinuous in the low-temperature region indicated by blue cross points. With increasing 
interband coupling, the tendency that the LLL modes of the order parameter describe the 
vortex lattice in high helds is enhanced. In the region where IA3I is on the order of A2, 
however, the hgure (b) shows that the Hc2 transition remains to be the second order even 
in a low-temperature limit [see the red circles in (b)]. This low-temperature second-order 
transition is diherent from the ordinary one occurring between the familiar vortex lattice and 
the normal phase in that the ordered phase below Hc2 is a vortex lattice with a PPB-induced 
modulation parallel to the held H || c. The fact that, in the present quasi-two-dimensional 
system with the simple cylindrical Fermi surfaces, a spatial modulation parallel to the c-axis 
tends to be created is surprising from the viewpoint based on the single-band model and 
should be regarded as one of the consequences of the multi-band electronic structure. Details 
of a possible modulated state in this case will be discussed in Sect. 4 . 

On the other hand, as IA3I becomes signihcantly large, the Hc2 transition again becomes 
discontinuous (mean-held hrst-order) in the low-temperature region indicated by the blue 
cross points, and no higher LL vortex lattice tends to appear even at lower temperatures. 
Furthermore, note that Fig. m c) belongs to the category with \U\ < 0 . Although we have 
not examined the details of the SC order parameter just below the discontinuous transition 
line Hc2{T), a spatial modulation parallel to the c-axis of the SC order parameter similar 
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to that of the LO state is expected to appear there. Note that the case with a sufficiently 
large interband coupling is qualitatively different from the single-band case. 


4. Summary and Discussion 


In this paper, we have attempted to understand the high-held superconductivity in quasi- 
two-dimensional two-band materials with a remarkable paramagnetic pair-breaking. Usually, 
the paramagnetic pair-breaking effect is expected to occur, even in two-band materials, in 


a situation with negligibly weak orbital pair-breaking 


23 | such as in the held conhgura- 


tion parallel to the conducting layers. In the present work, we have focused rather on the 
H II c case and have found new aspects on two-band superconductivity in high magnetic 
helds, which cannot be seen in the single-band superconductivity. In the single-band limit, 
the shape of the Hc2{T) curve and the character of the Hc2 transition depend only on the 
parameter a oc (y'/lvp proportional to the Maki parameter, which corresponds to the para¬ 
magnetic ehect relative to the orbital ehect, where g is the effective ^f-factor and v is the 
Fermi velocity. Thus, the mechanism of the pair-breaking depends only on the ratio of the 
paramagnetic effect relative to the orbital one. In two-band systems, however, a nonva¬ 
nishing interband coupling A3 induces new dependences of the Hc2 transition on material 
parameters. Not only the a-value on each band but also dependences on the parameters 
describing asymmetries between the two bands can change the nature of the high-held su¬ 
perconductivity. In other words, PPB can be regarded as acting independently of the orbital 
pair-breaking in two-band superconductors. For instance, a strong PPB is not equivalent to 
a weak orbital pair-breaking in these systems. 

To elucidate the parameter dependences of the phase diagram, we have changed the 
band asymmetries in the Fermi velocity (or the effective mass) and Zeeman energy (or the 
(/-factor) separately and have investigated their dependences on the interband coupling IA3I. 
The former band asymmetry primarily affects the orbital pair-breaking, while the latter 
measures a difference in PPB between the two bands. It is found that, when the band 
asymmetry in the Fermi velocity is dominant, the Hc2 transition at low temperatures tends 
to become a hrst order one, and that an additional PPB-induced spatial modulation in the 
vortex state below the hrst-order Hc2 transition can occur only in the direction parallel 
to the field and hence to the c-axis. On the other hand, when the band asymmetry and 
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thus the difference between the Maki parameters dehned on the two bands stems mainly 
from that in the Zeeman energy, we hnd that, for moderately large IA3I values, the Hc2 
transition at low temperatures remains to be continuous, while the PPB-induced spatial 
modulation just below Hr‘? at low temperatures can occur only along the held || c. As 
pointed out previously [l 9 |, these features correspond to the necessary conditions for the 
realization of the crisscrossing vortex lattice proposed in Ref. 20 . Therefore, a new high 
held SC phase can be expected to occur in quasi-two-dimensional two-band materials in 
H II c. On the other hand, when the Hc2 transition at low temperatures is of hrst order, 
the ordinary Larkin-Ovchinnikov state jo] modulating along the c-axis may occur just below 
Hc2{T). The present results might be relevant to the recent observation of an unexpected 


high held SC phase in FeSe 


18 | . To identify the observed high held phase in FeSe and, if any. 


the corresponding one in other iron-based superconductors with a theoretical one properly, 
further model calculations and studies on the huctuation ehect are needed. 

Appendix 


First, we explain how we use the HS transformation and how we obtain the GL functional 
in the form of Eq. (j3]). Next, the details of the derivation of the GL quadratic term ^2(2) 


and the quartic one (2(4) are explained, 
closely follows that in the previous work 


he treatment to derive (2(2) and (2(4) given below 


(lOj 


Hubbard-Stratonovich transformation for interaction term 


We show how we obtain the GL functional in the form of Eq. (| 3 ]). Here, we consider the 
case of | 17 | > 0 . By using the functional integral method, from the two-band Hamiltonian 
(Eq. (II])), we can write the partition function for grand-canonical ensemble as 


2:= 


( 20 ) 


where V’n(x, j, r) and the Grassmann numbers that represent fields of quasi¬ 

particles and r represents imaginary time. The noninteractive part of the action is given 
as 


^( 0 ) = [^ dr 

^ rt rr 


A 


+ ^n(-«V - eA(r)) - aln + J 


( 21 ) 
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where for simplicity, we define J as = Jl'ipnU + l)V’n(i) + U J + l)]/2- On 

the other hand, the interactive part is given as 

Sint = -fdTj2 J (22) 

^ j n,n' 

Now, to integrate formally qnasi-particles’ degrees of freedom, we introdnce anxiliary 
complex fields Ai(r,j, r) and A 2 (r,j, r) corresponding to the SC order parameter fields, 
and transform the interactive part in Eq. (1201) as 


,-Sint _ 


/3 


= / VA*VAexp \ - dTj2 dh 


EK(U-‘) An' 

n,n' 

- E + Ki'iii) 


(23) 


Hereafter, becanse we want to obtain the GL fnnctional that corresponds to the mean field 
theory, we neglect the qnantnm fluctuation or r-dependence of A„. As a result, we get the 
partition function as 


Z = / VA*VA exp 


PE f d^rEK{u-').^„A' 


VppVip [b,b]-5'[p,p,A*,A] ^24) 


nr riT n 


where S' = - dr J2j / d^r + A* 

Then, by using the linked-cluster theorem, we can transform Eq. fl24l) as 


Z = / VA*VA exp -13 


^(o)+E J K (r,j)(u A„,(r,j) 

j n,n' 

-T - l) |,(25) 

where r2(o) = — Tln[/D-^D-^exp(—S'(o)[t/’, V’])] fh^-f represents the thermodynamic potential 
in the normal phase and denotes the ensemble average of X with respect to the non¬ 

interactive part of the action with only the connected diagrams retained. Finally, expanding 
exp(—S") with respect to S' in Eq. (l25ll . we obtain the GL functional in the form of Eq. (|3]): 


fioLlA-.A] = fi(„)+i:/ d^,:EK{U-')nn'^’ 


n,n' 


- ^(s'y -'^(s'y +0(Ai6,A2®), 

91 \ /(O) 41 \ /(O') \ ^ ^ J ' 


(26) 


2! ' '(0) 4! /((I) 

where the minimum point of Ggl corresponds to an equilibrium state at the mean field level. 
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Calculation of Q 


( 2 ) 


First, iF( 2 )n(n, k) defined in Eq. ([6]) will be rewritten in a convenient form for nnmerical 
calculations. After integrating over and exponentiating the resulting denominator in 

terms of the integral 

1 r°° 

X ^ Jo 

where X > 0, we have 

A(2)„(n,fc) = 2TTTDn r dp ( 2^2 ^ LivnPP-u\ _^28) 


Here, a lower cutoff pcut, corresponding to a high energy cutoff, of the p-integral has been 
introduced. Then, to average exp{iVnP p-H) with respect to the solid angle (or p) in Eq. fl28l) . 
we introduce raising and lowering operators of the Landau level as 7 r± = A/f(na; ± iily)/^/2. 
By expanding the exponential, we obtain 


^ivr.pp-n\^ =exp 


VrJp^'\ ^ 

■4A^v 


Vn P \ -W-Af 
TT , 7r_ . 


2A 


(29) 


H 


From this equation, we see that the eigenstate of iF( 2 )n(n,/c) is in the n-th Landau level. 
Lastly, with respect to oji, a, and fc', we sum or average the functions expressed in Eq. (|28|) 
and then, with Eq. ([29]), we hnd 

/?(2l„(n,A^)7-„(r) = 2^TD„ r dp-^^^P^]J2Jpsmf] 

Jpcut smh(27rTp) \ 2/ 

“P {-Jx?) (if?) 

where Jbv(r) belongs to the A^-th Landau level, Jo(a() represents the zeroth-order Bessel 
function, and Ljv((c) represents the A^-th-order Laguerre polynomial. Now by substituting 
Eq. fl30l) in Eq. dSj), we get f2(2) in the form of Eq. ([7|). 


Calculation of H( 4 ) 

In this case, after performing the ^n(p) integral, the expression of iF( 4 )„({nj, fcj}) [Eq. 
flT^ ] consists of three denominators, and thus, the expression corresponding to Eq. fl28l) of 
A( 4 )„({nj, fcj}) is expressed as threefold integrals over the parameters pj {j = 1, 2, and 3). 
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After properly symmetrizing the expression w.r.t. Ylj and kj, A'( 4 )„({nj, ki}) takes the form 

1 ? _ 4vrTA)„ cos[24(pi + Pa + Ps)] t /^o t, 

-^(4)n({ni, A;^}) — / dpidp2dp^ . Jo (2 Jv-D) 

Jo sinh 27rT Pi + p 2 + P3 V / 


Alayer ^0 sinh [27rT(pi + P2 + Ps)] 

X (exp [iVnP ■ (pint + P 2 n 2 + PsH^)] )5, 


(31) 


where 


D({fc.,Pi}) = E 


i=l 


ki 

Pi sin I Y 


1 2 


+ 2 pip 3 Cos(^^^sin(^y^ 


+ 2 PiPjCos 

(i,i)=(l,2),(2,3) 

Pi\ . fkz 

sill 


sm — sm 


(32) 


and kij = ki — kj. 

Now, what we have to do is calculate the average of exp [iVnP ■ (piH^ + P2n2 + PsHg)] 
over the direction of p expressed in Eq. ([3T]). Below, we adopt the Landau gauge (A = Bxy) 
and assume that, on the transition line, the order parameter belongs to the LLL and has 
only component’ (see Sect. 2). Then, we can expand the order parameter as A„(r, k) = 
{SH)~^^‘^w^Pn{k',T, H) J2qUq{r)(j){q,k), where Sh = ^/TiXuLy with the system size Ly along 
the p-axis, Mq(r) = exp[—(a: + qXn'^)'^/{2X'jj) + iqy], and 0(g, k) is the expansion coefficient. 
We use the following identity: 

(|/^n,iP ~ {Xi/Xjj + qiXjJ + /3n,iY 


exp 


(ipiVnP ■ Il^Uq^{Ti) 


= exp 


+ ^q^yi 


, (33) 


4 2 

where /3n,i = PiVn{Px — iPy)/XH- We can show this by using the Campbell-Baker-Hausdorff 
formula. 

Finally, by substituting Eq. fIM]) in Eq. flTT|) . expanding the order parameter A„(rj,Pj) 
with Mg.(rj), and using Eq. fl55]l . we can explicitly conduct the integration in Eq. flTT]) with 
respect to x and p, and then obtain the hnal form of the quartic term 
1 

LLL ~ 


Q 


(4) 


Jit.ifcit Vi=l / 

cos[2/„(pi + P 2 + ps)] 


V^Ni^yerSn n,{qi},{ki} 

roo 

X (j)*{qi,ki)(j){q2,k2)(p*iq3,k3)(p{q4,k4)2TTTDn [ dpidp2dp3 
X Jo (2J'/D'j le 


P4=0 


sinh[27rT(pi + p 2 + ps)] 
(34) 


where qij = qi - qj, 


B{{qi,Pn,i}) = -- 


Y — {l^nP + fdn,2^ + 

1 


. 2=1 

4 


ECm -4Ec« 


\i=l 


2=1 


+ 7(^13 + 924 )Xh , 


(35) 
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Cn,i Pn,i 1) 3), cllicl (^n,i (^n,i 2, 4). 

Now, if we restrict the order parameter with modulation ±k along the c-axis, we obtain 
Eq. flT^ from Eq. (1541) . 
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Fig. 3. (Color online) Hc 2 {T) curves (black solid line with red circles or blue cross points) following 
from the same an values as in Fig. [2l Differing from Fig. [2l the band asymmetry is assumed here 
to arise primarily from the Zeeman energy: the parameters describing the breakdown of OnS are 
chosen as m\ = 2 m 2 and gi = 4 g (2 in this case. The intraband couplings are the same as those in 
Fig. [21 while only A 3 is changed as (a) 0.05, (b) 0.15, and (c) 0.4. For comparison, the N = 1 LL 
curve is also shown in each case as in Fig. [H and the transition curve, estimated roughly with no 
calculation, between the modulated high-held phase and the ordinary vortex lattice is drawn with 
the red dashed curve. In all cases, the Hc 2 {T) line is dehned in LLL, and the nature of the transition 
is continuous on the black solid curves. On the blue cross points in (a) and (c), the transition is 
discontinuous, while, on the red circles in (b), the Hc 2 transition becomes the second order to a 
vortex lattice with a modulation parallel to the held and c-axis. The pink-dashed curve expresses 
the temperature dependence of the magnitude of the modulation wavevector of the high-held SC 
phase on Hc 2 {T). 
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